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ABSTRACT 

In this paper we derive new expected scaling relations for clusters with giant radio 
halos in the framework of the re-acceleration scenario in a simplified, but physically 
motivated, form, namely: radio power (Pr) vs size of the radio emitting region (Rh), 
and Pr vs total cluster mass (Mr) contained in the emitting region and cluster velocity 
dispersion (an) in this region. 

We search for these correlations by analyzing the most recent radio and X-ray 
data available in the literature for a well known sample of clusters with giant radio 
halos. In particular we find a good correlation between Pr and Rh and a very tight 
"geometrical" scaling between Mjj and Rh- From these correlations Pr is also ex- 
pected to scale with Mh and <jh and this is confirmed by our analysis. We show that 
all the observed trends can be well reconciled with expectations in the case of a slight 
variation of the mean magnetic field strength in the radio halo volume with Mjj. A 
byproduct correlation between Rh and an is also found, and can be further tested by 
optical studies. In addition, we find that observationally Rh scales non-linear ly with 
the virial radius of the host cluster, and this immediately means that the fraction of 
the cluster volume which is radio emitting increases with cluster mass and thus that 
the non-thermal component in clusters is not self-similar. 

Key words: particle acceleration - turbulence - radiation mechanisms: non-thermal 
- galaxies: clusters: general - radio continuum: general - X-rays: general 



1 INTRODUCTION 

Radio halos are diffuse Mpc scales synchrotron radio sources 
observed at the center of a growing number (~ 20) of massive 
galaxy clusters (see e.g., Feretti 2005 for a review). Radio 
halos are always found in merging clusters (e.g., Buote 2001; 
Schuecker et al 2001) thus suggesting a link between the dy- 
namical status of clusters and the radio halos. Observations 
show that radio halos are rare; however present data suggest 
that their detection rate increases with increasing the X-ray 
luminosity of the host clusters and reaches 30-35% for galaxy 
clusters at z ^ 0.2 and with X-ray luminosity larger than 
10 45 erg/s (Giovannini, Tordi & Feretti 1999, GTF99). 

Two main possibilities have been so far investigated to 
explain the radio halos: i) the so-called re- acceleration mod- 
els, whereby relativistic electrons injected in the intra cluster 
medium (ICM) are re-energized in situ by various mecha- 
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nisms associated with the turbulence generated by massive 
merger events (e.g., Brunetti et al. 2001; Petrosian et al. 
2001); ii) the secondary electron models, whereby the rel- 
ativistic electrons are secondary products of the hadronic 
interactions of cosmic rays with the ICM (e.g., Dennison 
1980; Blasi & Colafrancesco 1999). 

Recently, calculations in the framework of the re- acceleration 
scenario have modelled the connection between radio ha- 
los and cosmological cluster mergers, and investigated the 
observed correlations between the synchrotron radio power 
and the X-ray properties of the hosting clusters (Cassano 
& Brunetti 2005, CB05; Cassano, Brunetti & Setti 2006, 
CBS06). Observed correlations relate the radio power at 1.4 
GHz (-Pi. 4) with the X-ray luminosity (Lx), temperature 
(T) and cluster mass (Liang 1999; Colafrancesco 1999; Fer- 
etti 2000,2003; Govoni et al. 2001a; Enfilin and Rottgering 
2002; CBS06); also a trend between the largest linear size of 
radio halos and the X-ray luminosities of the hosting clus- 
ters is found (Feretti 2000). In particular, CBS06 found a 
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correlation between Pi. 4 and the virial mass M v of the host- 
ing clusters, Pi. 4 oc M^' 9±0 ' 4 , by combining the Pi. 4 — Lx 
correlation derived from a sample of 17 giant radio halos 
with the M v — Lx correlation obtained for a large sample 
of galaxy cluster compiled by Reiprich & Boheringer (2002). 
However, this correlation, which has been discussed in the 
particle re-acceleration scenario by CBS06, relates quantities 
which pertain to very different spatial regions: the observed 
radio emission comes from a radial size Rh ~ 3 — 6 time 
smaller than the virial radius R v . 

In this paper we discuss expected scaling relations for radio 
halos in the framework of the re-acceleration scenario in its 
simplest form. Then, we derive a novel observed correlation 
between the radio power of radio halos and their extension 
and a tight "geometrical" correlation between the size of 
radio halos and the mass of the cluster within the emitting 
region. We also present additional correlations which are ex- 
pected on the basis of these two scalings. Finally we compare 
all these observed correlations with the model expectations. 

A ACDM (H = 70 km s" 1 Mpc" 1 , Q, m = 0.3, Ov = 
0.7) cosmology is adopted. 

2 PARTICLE ACCELERATION SCENARIO 

2.1 Main features and implications of the 
re-acceleration model 

The particle re-acceleration model is designed to explain the 
origin of the synchrotron radio emission diffused on scales 
larger than that of the cluster cores (giant radio halos), 
while the so called mini-halos and other smaller scale dif- 
fuse sources at the cluster center (e.g. core halo sources) 
might have a different origin (e.g., Gitti, Brunetti, Setti 
2002; Pfrommer & Ensslin 2004, and ref. therein). 

In the conventional particle re-acceleration scenario the 
lower energy electrons (7 ~ 100 — 300), relic of past activities 
in the clusters, are re-energized due to resonant and/or non- 
resonant interactions with the turbulence developed dur- 
ing cluster-cluster mergers. Turbulence and shear flows are 
expected to amplify the magnetic field in galaxy clusters 
(e.g., Dolag et al. 2002, 2005; Briiggen et al. 2005) however, 
the decay time-scale of the magnetic field is expected to be 
larger than several Gyr (e.g., Subramanian, Shukurov and 
Haugen 2006) and thus the particle re- acceleration process 
can be thought as occurring in a stationary magnetic field 
amplified during the previous merging history of the cluster. 

The basic features of this model can be briefly summa- 
rized as follows: 

a) The average synchrotron spectrum of radio halos is 
curved and can be approximated by a relatively steep quasi- 
power law which further steepens at higher frequencies up 
to a cut-off frequency. 

The curved, cut-off spectrum is a unique feature of the 
re-acceleration model, which well represents the typical ob- 
served radio halo spectrum, due to the existence of a max- 
imum energy of the radiating electrons (at j max < 10 5 ) 
determined by the balance between the energy gains (re- 
acceleration processes) and synchrotron and inverse Comp- 
ton losses (e.g., Brunetti et al. 2001, 2004; Ohno et al 2002; 
Kuo et al. 2003). Accordingly, the detection of a radio halo 
critically depends on cut-off frequency which should be suf- 



ficiently larger than the observing frequency. As a conse- 
quence, there is a threshold in the efficiency which should be 
overcome by the re-acceleration processes in order to accel- 
erate the electrons at the energies necessary to produce radio 
emission at the observed frequency in the clusters' magnetic 
fields. In the merger-related scenario it is expected that only 
mergers between massive galaxy clusters may be able to gen- 
erate enough turbulence on large scales to power giant radio 
halos at GHz frequencies, thus not all clusters which show 
some merger activity are expected to possess a giant radio 
halo. In particular, CB05 show that the expected fraction of 
clusters with radio halos increases with cluster mass due to 
a more efficient particle re-acceleration process in more mas- 
sive galaxy clusters, and this is in line with the increase of 
the fraction of radio halos with cluster mass which is claimed 
from the analysis of present radio surveys (e.g., GTF99). 

b) In the re-acceleration model radio halos should be 
transient phenomena in dynamically disturbed clusters. The 
time scale of the radio halo phenomena comes from the com- 
bination of the time necessary for the cascading of the tur- 
bulence from cluster scales to the smaller scales relevant for 
particle acceleration, of the time-scale for dissipation of the 
turbulence and of the cluster-cluster crossing time. 

Present observations suggest that radio halos are prefer- 
entially found in dynamically disturbed systems (e.g., Buote 
2001; Govoni et al. 2004). Under the hypothesis that radio 
halos form in merging clusters in the hierarchical scenario, 
Kuo et al. (2004) found that the lifetime of these radio ha- 
los should be < 1 Gyr to not overproduce the observed 
occurrence of these sources. 

2.2 Predicted scalings for giant radio halos 

In this Section we derive scaling expectations for giant and 
powerful radio halos in the context of the re-acceleration 
scenario in its simplest form. 

The most important ingredient is the energy of the tur- 
bulence injected in the ICM. Numerical simulations of merg- 
ing clusters show that infalling sub-halos induce turbulence 
(e.g., Roettiger, Loken & Barns 1997; Ricker & Sarazin 2001; 
Tormen, Moscardini & Yoshida 2004). An estimate of the 
energy of merging-injected turbulence has been recently de- 
rived in CB05 by assuming that a fraction of the PdV work 
done by the infalling sub-halos is injected into compressible 
turbulence. They show that the turbulent energy is expected 
to roughly scale with the thermal energy of the ICM, a result 
in line with recent analysis of numerical simulations (Vazza 
et al. 2006, V06). 

Once injected this turbulence is damped by Transit- 
Time-Damping (TTD) resonance with thermal and relativis- 
tic particles (at a rate T t h and P re ;, respectively). Since 
the damping time is shorter than the other relevant time 
scales (dynamical and re-acceleration) the energy density 
of the turbulence reaches a stationary condition given by 
£t/(Tth + Trei)) where it is the turbulence injection rate 
(CB05). When re-acceleration starts, the bulk of the en- 
ergy density of compressible modes which is damped by the 
relativistic electrons goes into the re-energization of these 
electrons. On the other hand, after a few re-acceleration 
times, in a time-scale of the order of the typical age of ra- 
dio halos, electrons are boosted at high energies at which 
radiative losses are severe (oc E 2 ) and the effect of parti- 
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cle re-acceleration (oc E) is balanced by that of radiative 
losses. The electron spectrum gradually approaches a quasi- 
stationary condition and it can be assumed that the energy 
flux of the turbulent modes which goes into relativistic elec- 
trons is essentially re-radiated via synchrotron and inverse 
Compton mechanisms: 
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where i ayn and ii C are the synchrotron and IC emissivities 
(and T th » Y re i, CB05; Brunetti & Lazarian 2007). 
The ratio iic/isyn simply depends on (Bcmb/Bn) 2 , where 
Bcmb = 3.2(1 + z) 2 fiG is the equivalent magnetic field 
strength of the CMB (z, the redshift) and Bh the mean 
magnetic field strength in the radio halo volume, which can 
be parameterized as Bh oc M^ 1 with Mh the total cluster 
mass within Rh (the average radius of the radio emitting 
region) . 

Based on CB05, the injection rate of the turbulence in 
the radio halo volume can be estimated as i t oc ~p H v 2 /T cros , 
where ~p H is the mean density of the ICM in the radio 
halo volume, Vi is the cluster-cluster impact velocity, «| oc 
M v /R v , and T cros oc (Ry/M v ) ' 5 is the cluster-cluster cross- 
ing time (see CB05) and is constant by definition of virial 
mass in the cosmological hierarchical model (e.g., Borgani 
2006, for a review). In the case Rh is larger than the cluster 
core radius it is vj oc M v /R v oc Mh / Rh and an, the veloc- 
ity dispersion inside Rh, is oh = G Mh / Rh ~ o 2 (for the 
sake of clarity in Fig. [1] we report a comparison between oh 
and o v for our sample of clusters with radio halos). Thus we 
shall simply assume that the injection rate of turbulence in 
the radio halo volume is it oc ~p H o H . The term T re i/T t h is 
oc e re i/e t h xVT (Brunetti 2006, Brunetti & Lazarian 2007), 
where T is the temperature of the cluster gas, and e r ei/eth 
is the ratio between the energy densities in relativistic par- 
ticles and in the thermal plasma. Although this ratio might 
reasonably vary from cluster to cluster, we shall assume that 
it does not appreciably change in any systematic way with 
cluster mass (or temperature), at least if one restricts to the 
relatively narrow range in cluster mass spanned by clusters 
with giant radio halos (see also the results from numerical 
simulations for cosmic rays in Jubelgas et al. 2006). Then 
from Eq[T]the total emitted radio power is: 
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where we have taken \ / T oc oh and T(z, Mh ,bn) = [l + 
(3.2 (1 + zf/B H ) 2 ] ■ The expression T (Fig[2j) is constant 
in the asymptotic limit B H » B 2 mb or in the simple case 
in which the rms magnetic field in the radio halo region is 
independent of the cluster mass. For B 2 H « B 2 mb one has 
that T~ x oc M 2 ^ H , thus in the general case the expected 
scaling is steeper (slightly for Bh of the order of a few /xG) 
than that obtained by assuming a constant T . 

It is important to stress here that the expression in Eqf2] 
is a general theoretical trend which implies sim ple scaling 
relations. Indeed, by taking oh ~ \J GMh / Rh and under 
the assumption that the mass scales with Rh as Mh oc Rfj 
(see also Sect. 3.2), Eq[2] (with T ~ constant) entails the 
correlations: 
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Figure 2. Function T 1 , normalized to the T 1 value for a mean 
M H = M m = 3 ■ 10 14 Mq, as a function of M H , for b H = 0.5 and 
assuming different values of the magnetic field Bh corresponding 
to the mean mass B m = 0.5, 1, 3, 6 fiG, from top to bottom. 
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the effect of a non constant T is a steepening (although not 
substantial for ~ [iG fields) of these scalings. 



3 OBSERVED SCALING RELATIONS IN 
CLUSTERS WITH RADIO HALO 

Motivated by the theoretical expectations outlined in the 
previous Section, we have searched for the predicted scaling 
relations in the available data set for giant radio halos. Op- 
eratively, we will first discuss the case of the Pr — Rh scaling 
expected in Eq(3] which will allow us to address the tricky 
point of the measure of Rh in radio halos, and then we will 
show that a tight observational Rh — Mh scaling exists for 
radio halos. Then, we will discuss and verify the byproduct 
observational scalings between Pr — Mh and Pr — oh- 

We consider a sample of 15 clusters with known giant 
radio halos (Rh ^ 300 kpc) already analyzed in CBS06, 
with the exclusion of CL0016+16, due to the lack of good 
radio images to measure Rh, and of A754, due to very com- 
plex radio structure. References for 14 giant radio halos are 
given in CBS06, while for A2256 we use the more recent ra- 
dio data from Clarke & EnBlin (2006). In TabfJJwe report 
the relevant observed and derived quantities for our sample. 



3.1 Radio power versus sizes of radio halos 

A direct scaling between Pr — Rh for radio halos is not 
reported in the literature. We want to check the existence 
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Figure 1. a) ag versus virial velocity dispersion for the galaxy clusters in our sample together with the best fit power-law &h cx <rj ; 
b) Circular velocity profile, <r(< r) = (GM{< r)/r)° , normalized to the virial value from Navarro-Frenk- White (NFW; Navarro, Frenk 
and White 1997) models with c = R v /r a =3,4,5 for the dashed, dotted and solid lines, respectively. The rectangle indicates the region 
of the radio halos: Rh /Rv ~ 0.16 — 0.3. <Ju/cr v varies by a factor of less than 15% (for a fixed c) in our sample. 



Table 1. In Col.(l): Cluster name. Col. (2): cluster redshift. Col. (3): logarithm of the radio power at 1.4 GHz, Pi. 4, in unit of Watt/Hz. 
Col. (4): logarithm of the size of the radio halos, Rh , in unit of kpc h~ ■ Col. (5): logarithm of the total cluster mass inside Rh, Mjj, 
in unit of solar masses. The references for the cluster redshift and radio power are reported in CBS06, while for A2256 we use the more 
recent radio data from Clarke & EnBlin (2006). 



cluster's 




z 


Log(Pi. 4 ) 


Log(P H ) 


Log(Mjj) 


Log(cr H ) 


name 






[Watt/Hz] 


[kpc h 7 g] 


[M h" 1 ] 


[km 2 s~ 2 ] 


1E50657-558 





2994 


25.45 ± 0.03 


2.84 ±0 


04 


14 


83 ± 0.07 


6.63 ± 0.08 


A2163 





2030 


25.27 ±0.01 


3.01 ±0 


04 


15 


02 ±0.05 


6.65 ± 0.07 


A2744 





3080 


25.23 ± 0.04 


2.90 ±0 


06 


14 


76 ±0.10 


6.49 ± 0.11 


A2219 





2280 


25.09 ± 0.02 


2.84 ±0 


05 


14 


66 ± 0.08 


6.46 ± 0.09 


A1914 





1712 


24.72 ± 0.02 


2.77 ±0 


04 


11 


68 ±0.05 


6.54 ±0.06 


A665 





1816 


24.60 ± 0.04 


2.84 ±0 


04 


14 


57 ±0.09 


6.37 ±0.10 


A520 





2010 


24.59 ± 0.04 


2.61 ±0 


04 


14 


21 ±0.10 


6.24 ±0.11 


A2254 





1780 


24.47 ± 0.04 


2.61 ±0 


03 








A2256 





0581 


23.91 ± 0.08 


2.63 ±0 


04 


14 


17 ±0.09 


6.18 ± 0.11 


A773 





2170 


24.24 ± 0.04 


2.71 ±0 


03 


14 


43 ±0.05 


6.36 ±0.06 


A545 





1530 


24.17 ±0.02 


2.58 ±0 


03 


14 


08 ± 0.30 


6.13 ±0.30 


A2319 





0559 


24.05 ± 0.04 


2.63 ±0 


02 


14 


30 ± 0.03 


6.30 ±0.03 


A1300 





3071 


24.78 ± 0.04 


2.76 ±0 


14 


14 


54 ±0.17 


6.42 ± 0.22 


Coma (A1656) 





0231 


23.86 ±0.04 


2.53 ±0 


01 


11 


12 ±0.03 


6.22 ± 0.03 


A2255 





0808 


23.95 ± 0.02 


2.65 ±0 


03 


14 


16 ±0.07 


6.14 ±0.07 



of a Pr — Rh correlation by making use of directly mea- 
surable quantities, such as the power and the radius at 1.4 
GHz. In the present literature it is customary to use the 
Largest Linear Size (LLS), obtained from the Largest Angu- 
lar Size (LAS) measured on the radio images as the largest 
extension of the 2<r or 3<r contour level, as a measure of 
the radio emitting region {e.g., Giovannini & Feretti 2000; 
Kempner & Sarazin 2001). Since a fraction of radio halos in 
our sample is characterized by a non-spherical morphology, 
meaning a non-circular projection on the plane of the sky, 
an adequate measure of a radio halo's size can be obtained 
by modelling the emitting volume with a spherical region 



of radius R H =VRmin x R ma x, Rmin and R m ax being the 
minimum and maximum radius measured on the 3cr radio 
isophotes. In this way we have derived the Rh values for 
all 15 radio halos, as reported in Tab[T] by making use of 
the most recent radio maps available in literature. In Fig0] 
we report Pi. 4 versus Rh for our sample. We find a clear 
trend with Rh increasing with Pi. 4, i.e., the more extended 
radio halos are also the most powerful. The best-fit of this 
correlation is given by: 



log 
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= (4.18 ± 0.68) log 



Rh 
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Figure 3. Pi. 4 versus Rh- The fit has been performed using 
a power-law form in the log-log space and the best fit slope is 
reported in the panel. 



(0.26 ±0.07) 



(6) 



A Spearman test yields a correlation coefficient of ~ 0.84 
and a s = 0.00011 significance, indicative of a relatively 
strong correlation. 



5f 
O 




2.4 2.6 2.8 3 3.2 

log(R H )[kpc] 

Figure 4. Ratios between the average surface brightness of each 
radio halo and the corresponding la noise level from the radio 
maps. The five most regular radio halos are earmarked. 



3.1.1 Uncertainties in the measure of the size of radio 
halos 

The dispersion of the P\a — Rh correlation is relatively large, 
a factor of ~ 2 in Rh , and this may be due to the errors as- 
sociated with the measure of Rh - Indeed, radio halos are low 
brightness diffuse radio sources which fade away gradually, 
until they are lost below the noise level of a given obser- 
vation. Thus, the measure of a physical size is not obvious 
and, in any case, it needs to be explored with great care. 
However, what is important here is not so much the precise 
measure of Rh for each radio halo, but rather the avoidance 
of selection effects which might force a correlation. 

In principle the sensitivity in the different maps may 
play a role because the most powerful radio halos are also 
the most bright ones (Feretti 2005), and thus they might 
appear more extended then the less powerful radio halos in 
the radio maps. To check if this effect is present, in Fig(4] 
we plot the ratios between the average surface brightness of 
each radio halo in our sample and the rms of each map used 
to get Rh- It is clear that there is some scattering in the 
distribution which would yield a corresponding dispersion 
in the accuracy of Rh, however, and most importantly, the 
ratios are randomly scattered, and there is no trend with 
Rh, i.e., fainter radio halos are usually imaged with a higher 
sensitivity and thus the Pi. 4 — Rh correlation cannot be 
forced by the maps used to derive Rh- 

An additional effort in assessing the reliability of Rh 
(and of the Pi. 4 — Rh correlation) would be to measure the 
radial brightness profile of regular radio halos which are not 
severely affected by powerful and extended radio sources. In 
our sample it is feasible to obtain accurate radial profiles 
from available data for the following radio halos: A2163, 



A2255, A2744, A545 and A2319. We take the data at 1.4 
GHz (Feretti et al. 2001, Govoni et al. 2005, Govoni et al. 
2001a, Bacchi et al. 2003, Feretti et al. 1997, respectively), 
and use the software package SYNAGE++ (Murgia 2001) 
to extract the radial brightness profiles, after subtraction of 
the embedded radio sources. 

In Fig[5]we report the integrated brightness profiles of 
these radio halos. It is seen that the profiles flatten with 
distance from the respective clusters centres, indicating that 
basically all the extended radio emission is caught and that 
it is possible to extract an accurate physical size. In Fig(6]we 
report for these 5 radio halos the comparison between Rh, 
estimated directly from 3a radio isophotes (see the above 
definition), and Pss and R75, i.e., the radii respectively con- 
taining the 85% and 75% of the flux of the radio halos. 
We apply the same procedure also to the case of the Coma 
cluster at 330 MHz for which a brightness profile and radio 
map were already presented in the literature (Govoni et al 
2001b). For Coma at 330 MHz we find R H ~ 520 h^ kpc 
and Rss ~ 610 h^Q kpc, which set Coma in a configuration 
similar to that of the other clusters in Fig(S] 

The linear, almost one-to-one correlation between Rh 
and i?85 and the relatively small dispersion, consistent with 
the uncertainties in the profiles due to source subtraction, 
prove that our definition of Rh is a simple but representative 
estimate of the physical size of radio halos. 

We note that the sensitivities of the radio maps, the 
physical sizes i?85 and powers Pi. 4 of the 5 regular halos 
are representatives of the values encompassed by the full 
radio halo sample. Moreover, for these 5 radio halos alone 
we find Pi. 4 oc Pg 5 25±0 ' 63 , fully consistent with the P1.4 — Rh 
correlation obtained for the total sample. 
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Figure 5. Integrated radial brightness profiles of the cluster radio halos in A2744, A2319, A545, A2163 and A2255 (from the top left to 
the bottom right corner). The errors in the profiles (which are in the range 5-10%) include the uncertainties in the sources subtraction 
and the statistical errors (note that in this integral presentation the errors are not independent). 




2.4 2.6 2.B 3 3.2 
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Figure 6. Radius enclosing the 85% (filled circles) and the 75% 
(open circle) of the total radio flux at 1.4 GHz obtained by the 
profiles (Fig[5} versus Rh estimated directly from the radio maps 
at 1.4 GHz. 



3.1.2 Possible biases in the selection of the sample 

Since the Pi. 4 — Rh correlation is the driving correlation, 
one has to check whether this correlation may not be forced 
by observational biases due to the selection of the radio 
halo population itself. Indeed the great majority of these ra- 
dio halos have been discovered by follow-ups of candidates, 
mostly identified from the NVSS which is surface brightness- 
limited for resolved source^ and this may introduce biases 
in the selected sample. 

The upper bound of the correlation is likely to be solid: 
objects as powerful as those at the upper end of the corre- 
lation (log Pi. 4 ^ 25) but with small Rh (similar to that of 
radio halos in the lower end of the correlation) should ap- 
pear in the NVSS up to the largest redshifts of the sample, 
since, even at 2 ~ 0.3, they should be JSlO times brighter 
than the low power radio halos in the correlation and ex- 
tended (~ 2.5'). As a matter of fact A545 (z=0.15) and 
A520 (z=0.2), which are among the smaller radio halos in 
our sample, are already detected in the NVSS up to a red- 
shift 0.2 and there is no reason why objects with similar 

1 The rms brightness fluctuations in the NVSS are 0.45 
mjy/beam (beam=45x45 arcsec, Condon et al. 1998); the NVSS 
is sensitive to radio sources with size < 10' — 15' (appropriate for 
radio halos at z > 0.05; GTF99). 
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Figure 7. Distribution of the slopes, S, of the Pi. 4 — Rh correla- 
tion obtained with our Monte Carlo procedure (with 400 trails). 
The arrow indicate the value of the observed best- fit slope ~ 4.18. 



Figure 8. Mh versus Rh for giant radio halos. The best-fit 
power-law and the value of the slope are also reported in the 
panel. 



extension, but ~ 8 — 10 times brighter than A545 and A520, 
should not have been detected at z ^ 0.3. 

The lower bound of the correlation deserves much care 
since the brightness limit of the NVSS may play some role. 
It is clear that present surveys may significantly affect the 
selection of the faint end of the radio halo population. How- 
ever, Feretti (2005) and Clarke (2005), have already con- 
cluded that the typical brightness of the powerful and giant 
radio halos are well above the detection limit. 

In any case, a brightness limit should drive a Pi. 4 oc R 2 H 
correlation, much flatter then the observed one. In order to 
provide a further compelling argument against observational 
biases, we have run Monte Carlo simulations. To this end 
we have randomly extracted brightness values of hypothet- 
ical radio halos within a factor of ~ 5 interval (consistently 
with the range spanned in our sample) above a given mini- 
mum brightness and each time randomly assigned Rh and z 
among the observed values. In Fig[7]we report the distribu- 
tion of the Pi. 4 — Rh slopes obtained with our Monte Carlo 
procedure and note that this distribution is peaked around 
~ 2.5 with a dispersion of ±0.4 (this is somewhat steeper 
than the expected Pi. 4 oc R 2 H due to the well known red- 
shift effect, however small given the small redshift range of 
our sample). The values of the slopes from the Monte Carlo 
procedure are far from the observed value (FigO and a sta- 
tistical test allows us to conclude that the probability that 
the observed Pia — Rh correlation is forced by observational 
biases is < 0.05%. 



3.2 Geometrical Mh Rh scaling for radio halos 

The existence of a possible tight scaling between the size of 
radio halos and the cluster mass within the emitting region is 
not reported in the literature. Yet an observational Mh—Rh 
scaling may be important to relate virial quantities 0% = 
GM V /R V (« GMh/Rh = fffr) with quantities (R H and 



Mh) which refer to the emitting region, and to test simple 
model expectations (Sect. 2.2). 

At this stage of the paper, the main difficulty concerns 
the measure of the cluster mass inside a volume of size Rh- 
Here the only possibility is to use the X-ray mass deter- 
mination based on the assumption of hydrostatic equilib- 
rium. Nevertheless, radio halo clusters are not well relaxed 
systems and thus the assumption of hydrostatic equilibrium 
and spherical symmetry may introduce sizeable errors in the 
mass determination. Several numerical simulation studies, 
which have been undertaken in order to determine whether 
the above assumptions introduce significant uncertainties 
in the mass estimates, indicate that in the case of merg- 
ing clusters the hydrostatic equilibrium method might lead 
to errors up to 40% of the true mass, which can be either 
overestimated or underestimated (e.g., Evrard et al. 1996; 
Rottiger et al. 1996; Schindler 1996; Rasia et al. 2006). This 
would cause an unavoidable scattering in the determination 
of the mass in our sample, although there are indications 
that a better agreement between the gravitational lensing, 
X-ray and optically determined cluster masses is achieved on 
scales larger than the X-ray core radii (e.g., Wu 1994; Allen 
1998; Wu et al. 1998), which is the case under consideration 
(Rh > r c ). 

However, what is important here is that the mass deter- 
mination does not introduce systematic errors which depend 
on the mass itself and which may thus affect the real trend 
of the Pi. 4 — Mh correlation. We thus compute the total 
gravitational cluster mass within the radius Rh as: 



M H = M tot (< Rh) 



3K B TR 3 H p 
firripG 



(7) 



where r c is the core radius, T the isothermal gas tempera- 
ture and f3 the ratio between the kinetic energy of the dark 
matter and that of the gas (/3-model; e.g., Sarazin 1986). 
We have excluded from our analysis A2254 for which no in- 
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Figure 9. Pi. 4 versus Mh for giant radio halos. The best-fit 
power-law and its slope are also reported in the panel. 
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Figure 10. Square of the velocity dispersion inside Rh versus 
Rh . The best-fit power-law and its slope are also reported in the 
panel. 



formation on the /3-model is available. For the remaining 
14 clusters references are given in CBS06. From Eq[7] one 
has that Mh oc Rh for Rh » r c and Mh oc R h for 
Rh « r c . In Fig[8]we plot Rh versus Mh for our sample: 
we find Mh oc R H 17±0 ' 19 , which falls in between the above 
asymptotic expectations. 




Log(R v [kpc]) 



Figure 11. Rh versus virial radius, of hosting clusters esti- 
mated from the Lx — M v correlation (see CBS06). In the panel is 
also reported the best-fit correlation. 

3.3 Radio power versus mass and velocity 
dispersion 

In principle, the two correlations discussed so far for giant 
radio halos, Pi. 4 — Rh and Mh — Rh, imply the existence of 
correlations between Pi. 4 — Mh and Pi. 4 — an- In particu- 
lar Pi. 4 is expected to roughly scale as M H ^~ 2 . In Fig[9]we 
report Pi. 4 versus Mh for our sample together with the best- 
fit: Pi. 4 oc AP^ 99 * ' 22 , which is indeed in line with the above 
expectation. A Spearman test of this correlation yields a cor- 
relation coefficient of ~ 0.91 and s = 7.3 • 10 -6 significance, 
indicative of a very strong correlation. 

Pi. 4 is expected to scale with <jh and we found for our 
sample a best-fit correlation: Pi. 4 oc (o"fj ) 4,64±1 ' or ; a Spear- 
man test yields a correlation coefficient of ~ 0.89 and to 
s — 2 • 10~ 5 significance, indicative of a very strong correla- 
tion. 

Finally, as a by-product of all the derived scalings, it 
is worth noticing that also a trend between Rh — oh is ex- 
pected f Fig llOf) . This finding might also be tested by obser- 
vations in the optical domain which can directly constrain 
the velocity dispersion. 



4 IMPLICATIONS OF THE DERIVED 
SCALINGS 

Given that the larger radio halos are also the most pow- 
erful ones and are hosted in the most massive clusters, we 
expect that the size of a giant radio halo should scale with 
the size of the hosting cluster. We estimate for each cluster 
of our sample the virial radius (R v ) by combining the virial 
mass-X-ray correlation (M„ — Lx', CBS06) and the virial 
radius- virial mass relation (e.g. Kitayama & Suto 1996). 
This method allows to reduce the effect of scattering due to 
the uncertainties in the mass measurements (and thus in the 
R v ) of merging galaxy clusters (see discussion in CBS06). In 
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Fig llll we plot Rh versus R v for our sample. The best fit 
gives Rh oc P^' 63±0 ' 50 , i.e. a pronounced non-linear increase 
of the size of the radio emitting region with the virial radius. 
A Spearman test yields a correlation coefficient of ~ 0.74 
and s = 0.0023 significance, indicative of a relatively strong 
correlation, albeit less strong than the others correlations 
found in this paper. 

Given that massive clusters are almost self similar (e.g. 
Rosati et al. 2002) one might have expected that Rh scales 
with R v and that the radial profiles of the radio emission 
are self-similar. On the contrary, our results prove that self- 
similarity is broken in the case of the non-thermal cluster 
components. This property of radio halos was also noticed 
by Kempner & Sarazin (200f), which used a sample of ra- 
dio halos taken from Feretti (2000) and found evidence for 
a trend of the Largest Linear Size, LLS, with the X-ray 

1/2 

luminosity in the O.f-2.4 keV band, LLS oc LJ , while a 
flatter scaling, LLS oc R v oc L^ 6 is expected in the case of 
a self-similarity. Their results imply Rh oc R v ; if one takes 
Rh ~ LLS, this is substantially in line with our findings, ft 
is also worth noticing that X-ray-radio comparison studies 
of a few radio halos indicates that the profile of the radio 
emission is typically broader than that of the thermal emis- 
sion (e.g., Govoni et al. 200fb). The two ingredients which 
should be responsible for the break of the self-similarity are 
the distributions of relativistic electrons and magnetic fields. 
In MHD cosmological simulations (Dolag et al. 2002, 2005) 
it is found that the magnetic field strength in cluster cores 
increases non- linearly with cluster mass (temperature) . This 
implies that the radio emitting volume should increase with 
cluster mass because the magnetic field at a given distance 
from the centre increases with increasing the mass. A de- 
tailed analysis of the magnetic field profiles of massive clus- 
ters from MHD simulations could be of help in testing if 
the magnetic field is the principal cause of the break of the 
self-similarity. 



5 PARTICLE RE-ACCELERATION MODEL 
AND OBSERVED SCALINGS 

Although we have been guided by the analysis of Eq[2] to 
predict the existence of scaling relationships, the observed 
correlations derived in Sec. 3 are actually independent from 
the form of this equation. To test Eq[2]against the observed 
quantities of our sample of radio halos we make use of the 
monochromatic Pi. 4 instead of the unavailable bolometric 
Pr. This is possible because the typical spectral shape of 
radio halos is a r ~ f.l — f.2 (P(^) oc v~ ar ) and thus the 
K-correction is not important (CBS06). 

In Fig lf 21 we report Pi .4 versus Mho\. The best fit 
gives Pi. 4 oc (Mh ojf) 1 - 24 "- 19 . The observed scaling is 
slightly steeper, but still in line with the linear scaling ex- 
pected from Eq[2]for T constant (dashed line). As already 
discussed in Sec. 2.2 T is constant for B H » B\ mh or in 
the case in which the rms magnetic field in the radio halo re- 
gion is quite independent from the cluster mass (small bn), 
while formally a non-constant T always implies a steep- 
ening of the Pi. 4 — Mh o h scaling. Namely, in the case of 
~ /LtG magnetic fields, by combining Eq[2]with the observed 
Mh — Rh correlation (Sec. 3.2, Fig[8]), one has that the 




Log(M H a3 [M km 3 s- 3 ]) 

Figure 12. Pi. 4 versus Mh& h . The best-fits correlations (solid 
line) and the predicted scaling with T ~ constant (dashed line) 
are reported. 



best-fit in Fig lf 21 is fulfilled by the model expectations for 
0.05 < b H < 0.39. 

In principle the fit can be used to set constraints on the 
values of the theoretical parameters entering the normaliza- 
tion of Eq(2] (namely ecii/eth, and the fraction of the PdV 
work which goes into turbulence), but we will not pursue 
this any further here (see CB05 for a discussion). 

It is important to stress that not only the trend in 
Fig. 1121 but also the existence of the correlations found in 
Sec. 3 could have been predicted on the basis of the re- 
acceleration model (Sec. 2, Eqs. [3] [4] [5j under the very rea- 
sonable assumption that Mh oc R H - Indeed, if one uses the 
observed scaling Mh oc P^ 17±0 ' 19 to fix the parameter a, 
from Eq(2j and assuming the most simple case in which T 
is constant, one finds Pi. 4 oc R H 9 and Pi. 4 oc Mh S , which 
are actually consistent (within the dispersion) with the ob- 
served correlations (Sec. 3); as in the case of the trend in 
Fig. 1121 an even better fulfillment of all these correlations is 
obtained for a slightly non-constant T . 

A relevant point which derive from the comparison 
of model expectation and observed correlations (unless 
B H » Bl mb ) is that, at least under our simplified ap- 
proach (Sec. 2.2), Bh does not critically depend on cluster 
mass inside Rh and that radio halos might essentially select 
the regions of the cluster volume in which the magnetic field 
strength is above some minimum value (say ~ fj,G level). It 
is important to note that a roughly constant Bh with clus- 
ter mass does not contradict the scaling of B, averaged in 
a fixed volume, with cluster mass (or temperature) found 
in the MHD simulations (B within the cluster core radius, 
r c ~ 300 hj^kpc), and also found in CBS06 (B averaged 
within a fixed region of ~ 720 h^g kpc size) , because the 
magnetic field Bh is averaged over a volume of radius Rh 
that becomes substantially larger than the core radius with 
increasing the cluster mass (Pf//r c goes from ~ f.f to ~ 3 
with increasing cluster mass in our sample). 
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6 SUMMARY & CONCLUSIONS 

The particle re-acceleration model is a promising possibility 
to explain the origin and properties of the giant radio halos 
(e.g., Brunetti 2004; Blasi 2004; Hwang 2004; Feretti 2005, 
for recent reviews). 

• In its simplest form, as assumed here (Sect. 2), it predicts 
a very simple relationship (Eq[2]) between the total radio 
power Pr, the total mass Mh within the radio halo, the 
gas velocity dispersion an and the average magnetic field 
Bh- Under the assumption of a tight scaling between 
Mh and the size Rh, and that the gas is in gravitational 
equilibrium, Eq[5] naturally translates into simple scaling 
relations: Pr - Rh, Pr- M h , and P R - a H (Eqs. H H |S). 

Motivated by the above theoretical considerations, we have 
searched for the existence of this type of correlations by 
analyzing a sample of 15 galaxy clusters with giant radio 
halos. A most important point here is the measure of the 
size Rh, in itself a non-trivial matter, since the brightest 
radio halos may appear more extended in the radio maps 
and this might force artificial correlations with radio power. 
A careful analysis of published 15 GHz radio maps of the 
radio halos of our sample shows that this effect is not present 
(Sec. 3. 1.1). From the same data set we derive a meaningful 
estimate of the radius for each radio halos. We also show 
that our procedure leads to estimates fully consistent with 
the measurements from the brightness profiles worked out 
from the data for the five most regular radio halos; this 
consistency holds over the total range spanned by Rh in 
our sample (Sect. 3.1.1). 

• We obtain a good, new correlation (correlation coeffi- 
cient ~ 0.84) between the observed radio power at 1.4 
GHz and the measured size of the radio halos in the form 
Pia oc _R^ 18±0 - 68 (Sect. 3.1). In Sect.3.1.2 we discuss in detail 
several selection effects which might affect this correlation 
and conclude that it is unlikely that the observed correlation 
is driven by observational biases. 

• We address observationally also the presence of a tight 
scaling between Mh and Rh and this allows us to relate 
virial quantities to quantities in the emitting region. 

• The presence of the Pr — Rh and Mh — Rh correlations 
implies also other correlations. We derive relatively strong 
correlations (Sect. 3.3) in the form: P 1A oc M^ 99±0 - 22 
and Pi. 4 oc (o"h ) 4 ' 64±1 ' 07 , and, as a byproduct, also 

2 „ r>0.90±0.25 

a h oc R H 

A correlation between the size Rh and the cluster virial 
radius, P„, is qualitatively expected in the framework of the 
particle re-acceleration model. 

• In Sec. 4 we compare Rh vs. R v for our sample of clusters 
with giant radio halos, obtaining the non-linear trend 
Rh oc i? 2 ' 63±0,5 °, i.e., the fraction of the cluster volume 
that is radio emitting significantly increases with the cluster 
mass. This break of the self-similarity, in line with previous 
suggestions (e.g., Kempner & Sarazin 2001), points to the 
changing distributions of the magnetic fields and relativistic 
electrons with cluster mass and, as such, is potentially 
important in constraining the physical parameters entering 
the hierarchical formation scenario, such as the turbulence 
injection scale and the magnetic field strength and profile. 
Finally, we note that, by combining the Rh — Rv and 



Pia — Rh correlations, one easily derives Pi. 4 oc M%, 
which is consistent with previous findings (Pi. 4 oc M 2 ' 9±0,4 ; 
CBS06). 

• These observed correlations are well understood in the 
framework of the particle re-acceleration model. Indeed, we 
show that the theoretical expectation (Eq(2]| is consistent 
with the data (see Fig ll2|l . Assuming a simple constant 
form for T in Eq[2] and the observed Mh — Rh scaling, 
which is necessary to fix the model parameter a (Sect. 2), 
the model expectations (Eqs. |4j [3] [5]l naturally trans- 
lates into Pi. 4 oc R H 9 , Pi. 4 oc M^ 8 and P1.4 oc (a H f A 
correlations, all consistent (within the dispersion) with 
the observed correlations; an even better fulfillment of all 
these correlations is obtained for a slightly non-constant 
T , which corresponds to « fiG field in the radio halo 
region. Unless it is B H » B^ mb , from the comparison 
of model expectations and observations we conclude that 
Bh should not strongly depend on Mh, and thus in our 
simplified scenario (Sec. 2.2) radio halos essentially trace 
the regions of « fiG fields in galaxy clusters in which 
particle acceleration is powered by turbulence. 

To conclude, the particle re-acceleration model, closely 
linked to the development of the turbulence in the hierar- 
chical formation scenario, appears to provide a viable and 
basic physical interpretation for all the correlations obtained 
so far with the available data for giant radio halos. Future 
deep radio surveys and upcoming data from LOFAR and 
LWA will be crucial to improve the statistics and to provide 
further constraints on the origin of radio halos. 



ACKNOWLEDGMENTS 

RC acknowledge the MPA in Garching for the hospital- 
ity during the preparation of this paper. We thank Matteo 
Murgia for the use of the SYNAGE++ program. We thank 
Luigina Feretti for providing the data for A2163, A545 and 
A2319, and Marco Bondi for useful discussions. This work 
is partially supported by MIUR and INAF under grants 
PRIN2004, PRIN2005 and PRIN-INAF2005. 

REFERENCES 

Allen S.W., 1998, MNRAS 296, 392 

Bacchi M., Feretti L., Giovannini G., Govoni F., 2003, A&A, 400, 
465 

Blasi P., 2004, JKAS 37, 483. 

Blasi P., Colafrancesco S., 1999, APh 12, 169 

Borgani S., 2006, Lectures for 2005 Guillermo Haro Summer 

School on Clusters, to appear in "Lecture notes in Physics" 

(Springer), astro-ph/0605575 
Briiggen, M., Ruszkowski, M., Simionescu, A., Hoeft, M., & Dalla 

Vecchia, C. 2005, ApJ, 631, L21 
Brunetti G., 2004, JKAS 37, 493 
Brunetti G., 2006, AN 327, 615 

Brunetti G., Setti G., Feretti L., Giovannini G, 2001, MNRAS 
320, 365 

Brunetti G., Blasi P., Cassano R., Gabici S., 2004, MNRAS 350, 
1174 

Brunetti G., Lazarian A., 2007, MNRAS in press, 
as tro-ph/070359T1 



New scaling relations in cluster radio halos and the re- acceleration model 11 



Buote D.A, 2001, ApJ 553, 15 

Cassano R. & Brunetti G., 2005, MNRAS 357, 1313 (CB05) 
Cassano R., Brunetti G., Setti G., 2006, MNRAS 

369,1577(CBS06) 
Clarke T.E., 2005, ASP Conference Series, Eds. by N.E. Kassim, 

M.R. Perez, W.Junor and P.A. Henning, Vol 345, pag. 227 
Clarke T.E., Enfilin T.A., 2006, AJ 131, 2900 
Colafrancesco S., 1999, in "Diffuse thermal and relativistic plasma 

in galaxy clusters". Bohringer H., Feretti L., Schuecker P. 

Eds., Garching, Germany: MPE 271, p.269 
Condon J. J., Cotton W.D., Greisen E.W., Yin Q.F., Perley R.A., 

Taylor G.B., Broderick J. J., 1998, AJ 115, 1693 
Dennison B., 1980, ApJ 239L 

Dolag K., Bartelmann M., Lesch H, 2002, A&A 387, 383 
Dolag K., Grasso D., Springel V., Tkachev I.,2005,JCAP 1, 9 
Evrard A.E., Mctzler C.A., Navarro J.F., 1996, ApJ 469, 494 
Enfilin T.A., Rottgering H., 2002, A&A, 396, 83 
Feretti L., 2000, Invited review at IAU 199 'The Universe at Low- 
Radio Frequencies' in Pune, India, 1999 
Feretti L., 2003, in 'Matter and Energy in Clusters of Galaxis', 
ASP Conf. Series, vol.301, p.143, eds. S. Bowyer and C.-Y. 
Hwang. 

Feretti L., 2005, in 'X-Ray and Radio Connections', published 
electronically by NRAO, eds. L.O.Sjouwerman and K.K.Dyer 

Feretti L., Giovannini G., Bohringer H., 1997, New Astronomy 2, 
501 

Feretti L., Fusco-Femiano R., Giovannini G., Govoni F., 2001, 

A&A, 373, 106 
Giovannini G., Tordi M., Feretti L., 1999, NewA 4, 141 
Giovannini G., Feretti L., 2000, New Astronomy 5, 335 in 

'Merging Processes in Galaxy Cluster', vol.272, p. 197, eds. 

L. Feretti, I.M.Gioia, G. Giovannini 
Gitti M., Brunetti G., Setti G., 2002, A&A 386, 456 
Govoni F., Feretti L., Giovannini G., Boringer H., Reiprich T.H., 

Murgia M., 2001a, A&A, 376, 803 
Govoni F., Enfilin T.A., Feretti L., Giovannini G., 2001b, A&A 

369, 441 

Govoni F., Markevitch M., Vikhlinin A., VanSpeybroeck L., Fer- 
etti, L., Giovannini C, 2004, ApJ 605, 695 

Govoni F., Murgia M., Feretti L., Giovannini G., Dallacasa D., 
Taylor G. B., 2005, A&A, 430, L5 

Henry J. P., Finoguenov A.; Briel U.G., 2004, ApJ 615, 181 

Hwang C.-Y., 2004, JKAS 37, 461 

Jubelgas M., Springel V., Ensslin T.A., Pfrommer C, 2006, sub. to 

AfcA, [astro-ph/0603485| 
Kempner J.C., Sarazin C.L., 2001, ApJ 548, 639 
Kitayama T., Suto Y., 1996, ApJ.469, 480 
Kuo P.-H., Hwang C.-Y.; Ip, W.-H., 2003, ApJ 594, 732 
Kuo, P.-H.; Hwang, C.-Y.; Ip, W.-H., 2004, ApJ 604, 108 
Liang H., 1999, in "Diffuse thermal and relativistic plasma in 

galaxy clusters". Bohringer H., Feretti L., Schuecker P. Eds., 

Garching, Germany : MPE 271, p.33 
Markevitch M., Mazzotta P., Vikhlinin A., Burke D., Butt Y., 

David L., Donnelly H., Forman W. R., Harris D., Kim D.-W., 

and 2 coauthors, 2003, ApJ 586, 19 
Murgia M, 2001, Ph.D Thesis, Univerity of Bologna 
Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ 490, 493 
Ohno H. Takizawa M., Shibata S., 2002, ApJ 577, 658 
Petrosian V., 2001, ApJ 557, 560 
Pfrommer C, Enfilin T.A., 2004, A&A 413, 17 
Rasia E., Ettori S., Moscardini L., Mazzotta P., Borgani S., Dolag 

K., Tormen G., Cheng L.M., Diaferio A, 2006, MNRAS 369, 

2013 

Rickcr P.M., Sarazin C.L., 2001, ApJ 561, 621 
Reiprich T.H., Bohringer H., 2002, ApJ 567, 716 
Rosati P., Borgani S., Norman C. 2002, ARA&A 40, 539 
Roettiger K., Burns J.O., Loken C, 1996, ApJ 473, 651 
Roettiger K., Loken C; Burns J.O., 1997, ApJS 109, 307 



Sarazin C.L., 1986, Reviews of Modern Physics 58, 1 
Schindler S., 1996, A&A 305, 756 

Schuecker P., Bohringer H.; Reiprich T.H., Feretti L., 2001, A&A 
378, 408 

Subramanian, K., Shukurov, A., & Haugen, N. E. L. 2006, MN- 
RAS 366, 1437 

Tormen G.; Moscardini L.; Yoshida N., 2004, MNRAS 350, 1397 
Vazza F., Tormen G., Cassano R., Brunetti G., Dolag K., 2006, 

MNRAS 369L, 14 (V06) 
Wu X.-P., 1994, ApJ 436, L115 

Wu X.-P., Chiueh T., Fang L.-Z., Xue Y.-X., 1998, MNRAS 301, 
861 



